Critical state analysis of orthogonal flux interactions in pinned superconductors 
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We show that, based on the critical state model for flux-line pinning in hard superconductors, one 
can assess the magnetic moment relaxation induced by the oscillations of a perpendicular magnetic 
field. Our theory follows a recent proposal of using phenomenological 2D modeling for the description 
of crossed field dynamics in high-T c superconductors arXiv: cond-mat/0703330 . 

Stationary regimes with either saturation to metastable configurations, or complete decay to 
the thermodynamic equilibrium are obtained. The transition between both types of response is 
related to the disappearance of a flux free core within the sample. As a common feature, a step-like 
dependence in the time relaxation is predicted for both cases. The theory may be applied to long 
bars of arbitrary and non homogeneous cross section, under in-plane magnetic field processes. 
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I. INTRODUCTION 

As it was emphasized by Brandt and MikitikJ^ mag- 
netic relaxation in superconductors is not always related 
to thermally activated flux creep. Thus, one can predict 
the magnetic moment relaxation induced by perpendic- 
ular ac field oscillations, within the critical state frame- 
work. Recall that such phenomenological approactP de- 
scribes the zero temperature limit of magnetic flux dy- 
namics in superconductors with pinning. The basic idea 
is that a given magnetic source produces some distri- 
bution of persistent current densities (critical currents) 
within the sample, determined by the actual excitation 
process through Faraday's law. 

Essentially, the model in Ref|T] ascribes the magnetic 
moment relaxation to internal variations of the current 
density distribution, related to the shaking effect from 
the transverse ac field, ft is of note that, within such 
a theory, and depending on the actual conditions of 
the shaking process, one can either predict relaxation 
toward the true equilibrium state of the superconduc- 
tor, or eventual freezing into some metastable configu- 
ration. This feature seems to be tightly linked to the 
critical state physics. However, some recent research 
about crossed field effects in high-T c superconductors^ 
has raised doubts about the generality of the com- 
plete/incomplete relaxation issues. 

Dimensionality is a relevant property in the trans- 
verse flux phenomenon, because, at least, two compo- 
nent magnetic systems are involved. Then, as it was 
remarked by Brandt and Mikitik, the appropriate crit- 
ical state problem to be considered is, minimally, two 
dimensional. In their work for thin samples, taking ad- 
vantage of the smallncss of the thickness /width ratio a 
quasi-one-dimcnsional analysis was given. Namely, they 
split the problem into two coupled one-dimensional flux 
pinning statements. 

On the other hand, at this point, the reader may re- 
call that critical state models for one-dimensional geome- 



tries (infinite slabs and cylinders) incorporating trans- 
verse flux effects have been issued in different forms 
However, the latter studies refer to magnetic field vectors 
with mutually perpendicular components, but parallel to 
the surface of the sample. This mathematical simpli- 
fication implies an important physical assumption: the 
dynamics of transverse flux lines includes cutting effects 
between them. ClerrPl has justified that, still, the phe- 
nomenon may be treated within the critical state frame- 
work. This is done via the introduction of the so-called 
parallel (to B) critical current density J c u . In contrast to 
the most familiar perpendicular critical current density 
J c ±, which relates to the flux pinning effect P J c \\ relates 
to the maximum angle between flux lines before undergo- 
ing a cutting process, ft has been shown that a number 
of experimental observations can be understood within 
the flux pinning/cutting framework. These include the 
so-called butterfly and collapse effects in the transverse 
magnetic moment 
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FIG. 1: Sketch of the experiment simulated in this work. An 
infinite bar of cross section 2w x d is subjected to the magnetic 
field excitation shown in the plot. H y is initially increased to 
the value H y a and then H x is cycled with periodicity r ac and 
amplitude H x<a . 
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As a breakthrough in the above scenario, Vanderbe- 
mden et alW have emphasized that, with moderate nu- 
merical effort, two-dimensional theories can be used for 
describing transverse flux effects. This rules out the ne- 
cessity of considering either very thin samples or incorpo- 
rating cutting phenomena. The modeling is realized by 
analyzing the evolution of in-plane magnetic fields within 
the cross section of a superconducting bar. In fact, it is 
apparent that the mutually perpendicular flux lines pene- 
trate basically along (and not across) each other. Within 
this picture, one can analyze transverse flux effects only 
related to the flux pinning mechanism (the J c ± limita- 
tion). However, the former investigations were done in 
a flux creep framework (power law E oc (J/J c ) n ), and 
for quite specific initial conditions. As a consequence, 
the authors did always observe complete relaxation by 
transverse shaking. In this work, we have studied the 
true critical state limit (n — > oo in the above E — J 
law) for two-dimensional transverse flux problems under 
a wide range of conditions. By using this ansatz, it will be 
shown that cither complete or incomplete flux relaxation 
may be expected. 

The article is organized as follows. First, in Sec [IT] wc 
describe the in-plane magnetic flux configurations under 
study and justify the application of a generalized Bean's 
model for the investigation of such phenomena. A numer- 
ical method, based on the concept of mutual inductance 
between circuits will be introduced. In Sec jIII| a set of 
numerical results of the actual magnetic flux dynamics 
in a number of crossed field configurations is presented. 
Based on the critical current distributions, and also re- 
lated to the structure of the penetrating field lines, we 
show the importance of the partial penetration concept 
in this problem. In close analogy to the one dimensional 
critical state problems, increasing the applied magnetic 
field progressively leads to the disappearance of the flux 
free region. Sec |IV| is devoted to analyze the relaxation of 
the magnetic moment, induced by oscillations of trans- 
verse magnetic fields. A systematic study for different 
amplitudes of the polarizing dc field, as well as of the os- 
cillating transverse component is presented. Decay to the 
true equilibrium magnetization is clearly related to the 
disappearance of the flux free core. The generality of the 
above properties is shown by application of the model to 
samples with inhomogeneities in the cross section. Thus, 
in Sec|V]we illustrate the evolution of the magnetic flux 
structure for samples with non homogeneous critical cur- 
rent density J c . 



II. CRITICAL STATE MODEL 

The basic configuration considered in the present work 
is sketched in Fig{T] A superconducting bar occupies the 
region defined by x\ < w, \y\ < d/2, \z\ < oo. The 
external magnetic field will be applied along the x and 
y axes. Starting from a zero field configuration, H y will 
be ramped to a final amplitude H y ^ a . Then H x will be 



cycled between the values ±H x ^ a in a linear fashion. We 
note that, corresponding to the hypothesis that the lower 
critical field of the superconductor may be neglected as 
compared to H x ^ a and H y ^ a , the equality B = [IqH will 
be used in what follows. In other words, the equilibrium 
magnetization of the sample is approximated by zero. 

The starting point for the application of the criti- 
cal state theory to the above problem is our variational 
statement^ within the {^4, J }-formulationP Such a rep- 
resentation in terms of the vector potential and current 
density allows to include finite size effects in a simple 
manner, and has been well elaborated in previous work. 
The equivalence to the more standard differential equa- 
tion statements based on the Maxwell equations is rigor- 
ously justified in RefJTl Thus, the quasi-stationary evo- 
lution of magnetic processes in a hard superconductor 
may be obtained from the constrained minimization of 
the functional 



T = 



J n+ l(f) • J n+ i (x ') _ J n (x) ■ J n +1 (x') 



+ — [ (Aj.h+I - A e<n ) ■ J a+ l{x) 
MO Jfi ' 



(1) 



Here, J n stands for the current density at the time layer 
nSt, f2 represents the superconducting region, and A e 
means the externally applied vector potential. Minimiza- 
tion is performed iteratively in time with </ n +i (a?) the un- 
known function for each step. As regards the constraints, 
recall that, for the present work, the minimization of T 
under the condition || J\\ < J c will be equivalent to the 
flux pinning criterion | Jj_| < J c ±. Indeed, by symmetry 
arguments one can show that J strictly flows along the 
bar, i.e.: J — (0,0, J z ) and then ||J|| = |Jj_|. 

On the technical side, an important advantage of the 
above mentioned problem's symmetry is the simplifica- 
tion of the numerical algorithm to be used for solving the 
variational statement. Thus, one can argue that the cur- 
rent density streamlines may be replaced by a collection 
of infinite straight wires, and the variational statement 
can be easily discretized to an algebraic form. The prob- 
lem is transformed into the minimization of the quadratic 
function 



+M0 J]] -fi,n+l(^c,n+l _ ^o,n) • 



(2) 



Here, the set of unknowns {/,} represents the collection 
of current lines flowing across the section of the sample, 
and My is their mutual inductance matrix. Minimization 
will be made under the set of constraints — I c < U < I c . 
The parameter I c is the critical current for an elementary 
wire, related to the flux pinning current density limita- 
tion. The mutual inductance coefficients will be evalu- 
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ated from the expressions 



Ma = 



Mi. 



Mo 

2tt [Xi 



i) 2 + (y t -yj) 2 



(3) 



which one can be obtain for a collection of parallel 
straight wires of circular cross section (radius a)P 

In summary, the calculation method is as follows. 
First, we compute the matrix My for a given grid of 
elements, describing the problem (most of the exam- 
ples presented later correspond to the rectangular case 
of I x m — 46 x 92 elements). Second, we solve the min- 
imization of the function in Eq.(j2l. Recall that one has 
I x m variables ({h}) constrained by I x m conditions. 
This process is done iteratively in time (n — 0,1,2,.. .) 
and one introduces the desired external bias, through the 
position dependent vector potential. In the case of a uni- 
form field with components (H Xl H y ,0) one can use the 
z— component A e = H x y — H y x. 

Finally, we want to mention some additional details 
that may be of help from the practical point of view. 
On the one side, the further use of symmetry consider- 
ations related to the homogeneity of the applied mag- 
netic field is to be advised. Thus, if one considers the 
rectangular cross section in the inset of Figjl] the in- 
version property I(x,y) = —I(—x,—y) may be recalled. 
Only one half of the variables have to be accounted for, 
and the computation simplifies noticeably. On the other 
hand, an important advantage of the preferred {A, J }- 
formulation is that the physical quantities of interest are 
obtained by integration of the former. Then, discretiza- 
tion errors are smoothed. On the contrary, if one uses 
an {i/}-formulation, J is obtained by differentiation and 
this results in the magnification of errors. 
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FIG. 2: (Color online) Critical current distributions within 
the cross section of the superconducting bar, subjected to 
the process in Figjl] The applied magnetic field vector is 
indicated by the central arrow. (+) indicates outward current 
flow (J c ), while (-) indicates inward current flow (— J c ). The 
central core is the subcritical region J = 0. The normalized 



values h x 



0.4 have been used for this plot. 



±J C ,0. The actual structure of critical/subcritical zones 
within the bar's cross section strongly depends on the 
{H x (t), H y (t)} process. Below, we describe some sys- 
tematic properties that have been observed. 



III. CROSSED FIELD DYNAMICS 

Below, we will present the calculated current density 
distributions and the corresponding magnetic field con- 
figurations induced by external excitations in the form 
illustrated in FigJT] For further analysis, we have per- 
formed calculations for a variety of values of the parame- 
ters H x a and H y ^ a . Both have been allowed to range well 
above and below the characteristic field H p , which deter- 
mines the full penetration regime for H along the x axis. 
Within the geometry considered in this work (w — d) one 
hasP H p ~ 0.85wJ c /2. 

Hereafter, dimensionless units for the magnetic field 
will be used, in terms of H p , i.e.: h = H/H p . 

A. Critical currents 

As one expects from theoretical considerations P the 
numerical solution of Eq.Q produces a current den- 
sity distribution with the unknowns Ii taking the values 



1. Partial penetration 

In Fig|2]we depict the critical current distributions in 
the initial part of the process illustrated in Figjl] To 
the upper left one can find the typical penetration pro- 
file for a uniform external field along one of the symme- 
try axes of the sample (Y-axis). In this case, this has 
been achieved when the condition h v a = 0.4 is met. The 
subsequent pictures correspond to a number of interme- 
diate steps within the first cycle for the h x component. 
This has been also ramped to the amplitude h x a = 0.4. 
The actual orientation of the applied field is indicated by 
the central vector, that is proportional to (h x , h y ) within 
each picture. Notice the similarity between the current 
free cores within this study and the analytical results in 
Ref JTOl for thin samples in oblique fields. 

Two facts are noticeable in the evolution of the current 
density penetration profiles: (i) the introduction of h, x in- 
duces upper and lower sheets of current, devoted to the 
shielding of the magnetic field variations, (ii) within the 
chosen range of parameters, one gets a central current 
free core. The core is initially distorted by the applica- 
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FIG. 3: (Color online) Same as Fig(2] but using the values 
hx, a = 1.6 and h y ,a = 0.4. Only the stationary regime after 
several cycles is shown. 



tion of h x , but rapidly acquires a stationary behavior. 
This property is apparent in the lower shots of Fig|2] 
where we display the evolution within the first half of 
the second cycle in h x . Recall that one obtains a dynam- 
ical structure for the ±J C regions, around a stable central 
core with I c = 0. 



2. Full penetration 

As one could infer from the knowledge of one- 
dimensional critical state problems, the subcritical region 
collapses to zero volume when the applied magnetic field 
reaches a certain intensity. In the two-dimensional case 
this is also true, but one can obtain a full penetration 
regime by combining {h xai h y ^ a } and the actual excita- 
tion process in a number of ways. In Fig|3]we show the 
current density distributions obtained for the stationary 
regime in the case h x>a = 1.6, h y a = 0.4. Recall that, 
together with the disappearance of the core, any trace 
of the initial critical state, related to the application of 
h y a is lost. The stationary critical current distribution 
is only related to the oscillating h x component. 



B. Magnetic flux configurations 

The resultant magnetic field profiles (applied plus in- 
duced) corresponding to the previous current density dis- 
tributions are shown in Figs|4j|6] The magnetic field 
lines have been obtained as the contour plots of the 
function A z (x,y), i.e.: they are given by the condition 
A z = constant. This is a straightforward property in 
the 2D geometry of this work. On the other hand, the 
expression 



A 1 



\/{x- Xi) 2 + (y- yi) 2 /a 



(4) 



1. Partial penetration 

Figj4] displays the magnetic field structure correspond- 
ing to the transient period in which the central core is es- 
tablished (upper rows). Recall that the magnetic field is 
excluded from the current free region (compare to Fig {2]) 
and also the trend of forming a stationary profile, sur- 
rounded by regions in which significant variations of the 
flux structure take place. Thus, in the lower rows, one 
can already observe that the inner core is unperturbed, 
while the applied field oscillates in between ±h x a . 

Fig|4] has been obtained by integration of the current 
density profiles presented in the previous section, and 
thus corresponds to the crossed field amplitudes h y , a — 
0.4 and h xa = 0.4. 



2. The flux free core 

In order to have a more systematic information about 
the partial penetration regime in 2D problems with 
crossed magnetic fields, we have studied the formation 
of the flux free core in different conditions. Thus, Fig(5] 
displays the stationary core that is obtained when one 
uses an oscillating field amplitude h x<a — 0.4, for initial 
ramps with either h y>a — 0.2 or h v>a — 0.4 respectively. 




has been used in order to evaluate the vector potential 
contributed by an infinite wire of radius a, and carrying 
a current L. 



FIG. 4: (Color online) Magnetic flux lines for the process in 
Fig[2] The flux free core is defined by the two lines with the 
lowest magnetic field modulus within our numerical resolu- 
tion. 



5 



The behavior for other values has been also studied, and 
the results (not shown for brevity) do not display signif- 
icant differences. 

For the sake of clarity, in the plots, we have only shown 
the lowest level magnetic field lines that define the core, 
and their evolution as h x is cycled within the stationary 
regime. Notice that the boundary of the core remains 
unchanged, while the flux lines adapt to the applied field 
as one moves away from the center of the sample. The 
only difference between the two situations is the actual 
size of the core. As one can deduce from the results in the 
following section, this is straightforwardly related to the 
saturation value m y acquired by the vertical magnetic 
moment. 



Full penetration 



As it was shown above (Sec III A 2 1, the current free 



core disappears when one uses the external conditions 
hy.a — 0.4 and h x ^ a = 1.6. Here we show the flux line 
structure corresponding to the oscillations of h x in such 
a range. As one can see in Figj6] magnetic flux has fully 
penetrated the sample. However, the over damped char- 
acter of flux motion in the critical state is evident in that 
plot. Although the flux lines close to the periphery of the 
sample tend to follow the external excitation, the inner 
regions always presents a delay. 



IV. MAGNETIC MOMENT RELAXATION 

The connection between the magnetic field profiles ob- 
tained in the previous section, and direct experimental 
observations is established below. As a figure of merit 
we have chosen the sample's magnetic moment (per unit 
length) both along the x and y axes. The definitions 



d.r 



dx 



d/2 



-d/2 
d/2 



dy [yJ z (x,y)} 



d/2 



dy [xJ z (x,y)} 



(5) 



will be used. These quantities are straightforwardly cal- 
culated from the current distributions {J^} obtained by 
the method in Sec|TT| In practice, m x and m y can be 
measured by flux detecting coils, conveniently oriented 
around the sample. 

Next, we discuss the influence of h x ^ a and h ya on m x 
and m y , so as to provide testing criteria for our theory, 
and also to allow comparison with the available literature 
on the subject. 



A. Influence of the oscillation amplitude (h x ^ a ) 

First, we concentrate on the effect of changing h x _ a 
with h Vt a fixed, i.e.: one polarizes the sample by means 





x,a 



FIG. 5: Evolution of the magnetic field lines that define the 
stationary flux free core along the ^h x ^ a cycle. The upper plot 



corresponds to the vertical amplitude h y 



0.2, while the 



lower one was obtained for h v 



0.4. Dashed, dot-dashed, 



dotted and continuous styles are used for the consecutive cen- 
tral line pair during the cycle. 



of a given value of the vertical field, and performs etc cy- 
cles of the horizontal component for different amplitudes. 
According to our critical state theory, the results of such 
experiment would be as shown in Fig(7] Recall that, on 
taking h y a — 0.4 and cycling h xa between the values 
±h XtU = ±0.1, ±0.2, ±0.4, ±0.8, ±1.2, ±1.6, m x and m y 
display a behavior which is clearly different if one probes 
cither above or below h x a = 0.8. 






FIG. 6: (Color online) Magnetic flux lines corresponding to 
the full penetration profiles depicted in Fig|3] As before, the 
central vector is proportional to the applied magnetic field. 
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FIG. 7: (Color online) Evolution of the magnetic moment components for several excitation processes in the form sketched in 
Figjl] m x and m y are obtained through Eq.|5| and normalized to the values of m y at the beginning of the ac cycle. Time (t) 
is given in units of to = T ac w/2d. h x is given in units of the penetration field H v . The different lines (some are labeled for 
clarity) correspond to h x<a — 0.1, 0.2, 0.4, 0.8, 1.2, 1.6. All have been obtained for h y<a = 0.4 



As relates to m Xl the behavior for low amplitudes 
is characterized by a non-saturated oscillation that es- 
sentially follows the applied field h x . When plotted 
against this quantity m x reaches a stationary Bean's like 
loop after an initial transient in which the loop does 
not close. This behavior was already reported in both 
experimental and theoreticaP works. In the latter, by 
contrast to the present case, the flux cutting mechanism 
was considered. The response to higher amplitudes of 
the oscillating field (h xa > 0.8) is characterized by a 
saturation of m x , which relates to a stationary m x (h x ) 
loop, in the typical form of ID problems, when the full 
penetration regime is reached. 

The properties of m y are of special mention, since we 
have found less common features. Firstly, we outline the 
step-like behavior in the plot m y (t). Note that the de- 
scent of this quantity, induced by the oscillation of h Xl 
presents a series of plateaus that are precisely triggered 
by the return points in ±h x a . This fact is also appar- 
ent within the m y (h x ) plot, and was mentioned before in 
Refill where a power law E — J model and the high field 
region (full penetration) were considered. Recall that 
here, we also obtain the steps for the low field region. 
Nevertheless, the early saturation of m y makes the effect 



less visible in such cases. 

We want to stress the high similarity of our plot m y (t) 
and Fig. 2 in Ref. 1 . Thus, the main feature of both graphs 
is the separation of complete and incomplete relaxation 
when the value of h x ^ a approaches one in units of H p . On 
the other side, recall that the step- like behavior in our 
plot would be smeared out over a time scale comparable 
to the one used by the authors in Ref[TJ In fact, if one 
takes the time unit as to = T ac w/2d, the approximation 
for very thin samples used by Brandt and Mikitik leads to 
to S> T ac /2, while in our case (u> = d) we have to = r QC /2. 

B. Influence of the dc field amplitude (h Vl a) 

Next, we analyze the relevance of the initial h y ramp 
on the subsequent complete or incomplete relaxation of 
m y during the h x cycles. The main results for this study 
are plotted in FigjS] There, we show the evolution of 
m y for the oscillation amplitudes h x , a — 0.2,0.4,0.8,1.6 
and three different values of the dc field reached in the 
initial stage of the process: h Vta — 0.2,0.4,0.8. We em- 
phasize two aspects: (i) on the one side, the higher h Via , 
the lower value of the relaxed component m V)00 when re- 
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FIG. 8: (Color online) Time dependence of the normalized 
magnetic moment m y at various amplitudes of the ac mag- 
netic field h Xl a = 0.2, 0.4, 0.8, 1.6. For each value of h x>a we 
compare the decay of m y at three values of h y . a : h Vj a = 0.2 
(continuous line), h y>a = 0.4 (dashed line), h y , a = 0.8 (dot- 
dashed) 

laxation is incomplete. However, the separation between 
the values of m 9 CO first increases, and then decreases, 
going to zero when the complete relaxation is reached; 
(ii) on the other side, crossings between the values of the 
normalized magnetization are observed during the initial 
transient evolution. The importance of these crossings 
increases as the value of h y ^ a does. 



V. EXTENSIONS OF THE THEORY 

As an advantage of the mathematical modeling pro- 
posed in this work, one can deal with fully arbitrary 
2D problems, including non-uniform magnetic sources 
and non homogeneous properties in the cross section of 
the sample. The {A, J }-formulation stated in Eqs.|l]) 
and ([2| allows to introduce irregular cross sections just 
by defining the positions of the appropriate elementary 
wires, and calculating the corresponding Mij matrix. 

In this section, we show the modifications introduced 
by a non homogeneous current carrying capacity. First, 
we calculate the evolution of the magnetic field profiles 
for a sample with a central hole, when subjected to the 
process in Fig[T] Second, we consider the influence of 
having a central region with J c noticeably above the cor- 
responding value at the periphery of the sample. 



A. Samples with holes 

A sample with a hole is straightforwardly treated 
within the formulation in Eq.Q, just by skipping the 
variables related to the elementary wires within the 
empty region. Results are shown in Fig|9] This plot 
has been obtained for the values of the applied field 
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FIG. 9: (Color online) Magnetic flux lines corresponding to 
the excitation process in Fig[l]for a superconducting bar with 
a hole (indicated by the dashed line). We show several steps 
at the transient process within the first cycle. 



h-x.a = h y ,a = 0.4. Several stages at the first ac cycle 
are shown, which already suggest the formation of a sta- 
tionary regime, as in the case of Fig|4] Nevertheless, in 
the present situation, the transverse field shaking process 
cannot induce penetration of current within the empty 
region. On the other side, the magnetic field lines pro- 
gressively enter. Then, as a topological effect, the current 
free and the flux free regions are no longer coincident, 
subsequent to the contact between the penetrating flux 
and the hole. 

Just for the sake of brevity, we do not show the sta- 
tionary oscillations which follow after the last stage in 
Figj9] (one can easily guess the result). As for the case 
of homogeneous samples (Figj4j) they are characterized 
by a frozen core and swinging flux lines in the peripheral 
region. 




FIG. 10: (Color online) Same as Fig[9]for a full superconduct- 
ing bar with inhomogeneous pinning properties. J c within the 
dashed region is 5 times larger than around. 
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B. Non homogeneous critical current 

Fig |10| shows the evolution of the flux penetration pro- 
files for a sample with a higher pinning force in the 
central region. To be specific, we have used the ratio 
Jc,in = §Jc,around in the simulation, i.e. the constraint 
for the corresponding elements is — 51 c < I{ < 5I C . In di- 
mcnsionless units, the applied field amplitudes are again 

hx,a — tly.a 0.4. 

As one could expect, both flux and current penetra- 
tion into the central region are hindered by the higher 
pinning force. Note that when the initial flux free core 
is distorted by the oscillation of h Xl its boundary slips 
along the separation between the two regions. On the 
other hand, the establishment of the stationary regime 
basically consists of a vertical critical state profile within 
the high J c region, and the surrounding swing of flux 
lines described in the previous examples. 



VI. CONCLUDING REMARKS 

Although the majority of theoretical and experimen- 
tal studies about flux pinning in superconductors have 
focused on the dynamics induced by single-component 
applied magnetic fields, increasing interest in multi- 
component s ituations is arising. Along this line, some 
recent papers^- 3 * 4 * 5 * 11 * 12 ! have considered the importance 
of transverse flux effects. Of particular interest is the 
magnetization decay induced by the oscillations of a per- 
pendicular magnetic field. In this work, we have tried to 
clarify some questions about the ability of phenomeno- 
logical flux pinning models for describing the experimen- 
tal observations. In particular, we have addressed the 
debate about the existence of complete and incomplete 
relaxation regimes (full decay or not into the equilibrium 
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FIG. 11: A flux line moves to the right, by pivoting around 
the positions marked (+) and (-), when the ac field undergoes 
a cycle (sequence of numbers 1 — 6). The flux line has been 
selected from the set of profiles obtained after the conditions 
in Fig[6] The picture is split in two panels for clarity. 



magnetization) in terms of the actual transverse shaking 
process. We have shown that, within the critical state 
model for flux pinning, the kind of relaxation is related to 
the existence of a flux free core within the sample. Thus, 
under certain conditions (small enough amplitudes of the 
applied field components) one finds a transient regime in 
which the central flux free core changes its shape until a 
stationary profile is reached. Afterwards, the transverse 
field shaking produces the swing of the flux lines around 
the frozen core structure, and the magnetic moment be- 
comes constant (but not zero). If the amplitude of the 
applied magnetic field increases enough, the flux free core 
shrinks to null volume, and the stationary magnetic mo- 
ment becomes zero. 

The existence of incomplete/complete relaxation can 
be observed in experimental works^ and is also predicted 
by the vortex-shaking model in Ref , 1 . Although there are 
some differences between the statement of that model 
and ours, the essential implications are easily reconciled. 
To start with, the model in Ref|I]is a good approximation 
for thin strips, which are nothing but a limiting case of 
the 2D problem. On the other hand, it is of mention that 
a fundamental ingredient of Ref[T] is the idea that flux 
lines walk towards the center of the sample (x = 0), piv- 
oting around the so-called swivel points. Switching from 
one swivel point to the next is done every half-cycle. In 
the present work, we have merely used the concept of a 
maximum pinning force (introduced by the phenomeno- 
logical parameter J c ± , through the restriction | Jj_ | < J c ± 
in the variational solution of Maxwell equations). Re- 
markably, the evolution of flux lines in our model follows 
the walking pattern introduced by Brandt and Mikitik. 
This is shown in FigfTTj Starting from the magnetic field 
profiles related to the process in Fig|6j we have plotted 
a sequence of pictures, corresponding to a specific flux 
line, along the first cycle. It is apparent that two swivel 
points are formed (labeled + and — ), which establish the 
shift from position 1 to 6, by two rotations. The switch 
from the pivot + to — is clearly established by the re- 
turn point —H x a (step 3 to 4). We want to note that 
the vertical position of the swivel points has changed. 

Another point to consider is the role of the aspect ra- 
tio parameter w/d. Thus, the main difference between 
the predicted magnetization decay, that is, the step- like 
structure in our FigsjT] and [8] disappears when the phys- 
ical time unit to = r ac w/2d is taken into account. The 
steps are smeared out by the time scale of the plot, and 
the two models coincide again. 

Finally, a comparison between our proposal with the 
method introduced in RefJ3] (E cx (J/J c ) n ) also seems 
to be of interest. In principle, the results of this article 
should be obtained as the limit of the former, when higher 
and higher values of n are used. This would allow to 
quantify the importance of flux creep effects, intrinsically 
included in phenomenological E — J models. 
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